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Abstract 



I present a unifying scheme for hadronic functions that comprises logarith- 
mic corrections due to gluon emission in perturbative QCD, as well as power- 
behaved corrections of nonperturbative origin. The latter are derived by de- 
manding that perturbatively resummed partonic observables should be ana- 
lytic in the whole Q 2 -plane if they are to be related to physical observables 
measured in experiments. I also show phenomenological consequences of this 
approach. The focus is on the electromagnetic pion form factor to illustrate 
both effects, Sudakov logarithms and power corrections in leading order of 
Aq CD /Q 2 . The same approach applied to the inclusive Drell-Yan cross sec- 
tion enables us to perform an absolutely normalized calculation of the leading 
power correction in ^ 2 Aq CD (b being the impact parameter), which after ex- 
ponentiation, gives rise to a nonperturbative Sudakov-type contribution that 
provides enhancement rather than suppression, hence partly counteracting 
the perturbative Sudakov suppression. 
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1. INTRODUCTION 



In recent years, effort in QCD has turned increasingly toward the problem of including 
resummation effects due to multiple soft gluon emission, both in perturbation theory, as 
well as in the nonperturbative regime. The first effect is related to Sudakov suppression [1], 
well-known from QED, whereas those in the nonperturbative regime manifest themselves as 
power-behaved corrections [2], which, after exponentiation, amount to a Sudakov-like form 
factor [3]. However, as it turns out [4] this contribution provides enhancement rather than 
suppression. The hope is that improving the perturbative and nonperturbative structure of 
the theory this way, it will be possible to get better agreement with the existing hadronic 
data in terms of both correct overall shape and also normalization. In these investigations 
the crucial organizing principle is QCD factorization, which provides a handle to separate the 
short-distance (hard) component of a reaction (controlled by the large mass scale in the pro- 
cess, Q) - that will be treated perturbatively - from its long-distance (soft) nonperturbative 
part, related to the nontrivial QCD vacuum structure (and field condensates). 

In processes which involve the emission of virtual gluon quanta of low momentum, one 
must resum their contributions to all orders of the strong coupling constant. This gives rise 
to exponentially suppressing factors in 6-space (where b is the impact parameter conjugate 
to the transverse momentum Q±) of the reaction amplitude (or cross-section) of the Sudakov 
type with exponents containing double and single logarithms of the large mass scale of the 
process [1]. However, because of the Landau singularity of the running coupling at transverse 
distances b oc 1/Aqcd , an essential singularity appears in the Sudakov factor. Thus, one has 
to consider power corrections of O (& 2 Aq CD ) , which, though negligible for small b relative 

to logarithmic corrections oc In (& 2 Aq CD ) , may become important for larger values of the 
impact parameter. 

In this talk, I will discuss a general methodology to treat (power) series in the running 
strong coupling in connection with gluon emission. To be more precise, I will address this 
issue in terms of two processes: one to which the OPE applies, viz. the pion electromagnetic 
form factor at leading perturbative order, and another, the Drell-Yan process, to which the 
OPE is not applicable. The first is a typical example of an exclusive process with registered 
intact hadrons in the initial and final states (for a recent review and references, see, e.g., [5]). 
Such processes provide a "window" to view the detailed structure of hadrons in terms of 
quarks and gluons at Fermi level (Hadron Femptoscopy). The Drell-Yan mechanism, on the 
other hand, has two identified hadrons in the initial state and a lepton pair (plus unspecified 
particles) in the final state, whose transverse momentum distribution is proportional to the 
large invariant mass of the materialized photon. 

The goal in the second case will be to obtain not only the usual resummed (Sudakov) 
expression (which comprises logarithmic corrections due to soft-gluon radiation), but also 
to include the leading power correction as well, specifying, in particular, its concomitant 
coefficient. This becomes possible within a theoretical scheme, which models the IR behavior 
of the running coupling by demanding analyticity of physical observables (in the complex 
Q 2 plane) as a whole - as opposed to imposing analyticity of individual powers, i.e., order by 
order in perturbation theory -, while preserving renormalization-group invariance (references 
and additional information can be found in the recent surveys [6,7] and D.V. Shirkov, these 
proceedings). The underlying idea behind our method [4], is to demand that if hadronic 
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observables, calculated at the partonic level, are to be compared with experimental data, 
they have to be analytic in the entire Q 2 plane. This "analytization" procedure encompasses 
Renormalization Group (RG) invariance (i.e., resummation of UV logarithms and correct 
UV asymptotics) and causality (which imposes a spectral representation). As we shall see 
below, analytization removes all unphysical singularities in the the IR region, rendering 
perturbatively calculated hadronic observables IR-renormalon free. 



2. ANALYTIC FACTORIZATION SCHEME (AFS) 
2.1 Perturbative Pion Form factor with Sudakov Corrections 



F* (Q 2 ) 



pi poo rl^Yi 

^ dxdy j ^ VT (y, b, P'; C u C 2 , C 4 ) T H (x, y, b, Q; C 3 , C 4 ) 



Let us conduct our investigation by considering the space-like electromagnetic pion's 
form factor in the transverse (impact) configuration space: 

- 1 , , d 2 h 

xv i »(x,b,P;C 1 ,c 2 ,c A ) + ..., (i) 

where the modified pion wave function is defined in terms of matrix elements, viz., 

V n (x,b,P,fi) = |' k±l< Vk ± e^- b ^(x,k ± ,P) 

r dz ixP + z 



2vr 



e~ tx ^ z (0 



T (g(0) 7 + 7 5 g (o, z-,b ± )) I *( p )) A+=0 (2) 



with P + = Qj\pl = P~', Q 2 = —{P' — P) 2 , whereas the dependence on the renormalization 
scale /i on the RHS of (2) enters through the normalization scale of the current operator, 
evaluated on the light cone, and the dependence on the effective quark mass has not been 
displayed explicitly. In (2), T H is the amplitude for a quark and an anti-quark to scatter 
via a series of hard-gluon exchanges with gluonic transverse momenta (alias inter-quark 
transverse distances) not neglected from the outset. In the above, the ellipsis indicates 
the non-factorizing soft part, as well as disregarded higher-order corrections. The scheme 
constants C, emerge from the truncation of the perturbative series and would be absent 
if one was able to derive all-order expressions in the coupling constant. The scale C\jb 
{C\ = C3) serves to separate perturbative from non-perturbative transverse distances (lower 
factorization scale of the Sudakov regime and transverse cutoff). The re-summation range in 
the Sudakov form factor is limited from above by the scale C 2 ^Q (upper factorization scale 
of the Sudakov regime and collinear cutoff). 1 The arbitrary constant C4 serves to define 
the renormalization scale Cif{x 1 y)Q = //r, which appears in the argument of the analytic 
running coupling af 1 [8] (choice of renormalization prescription): 

-an(l) ( g 2) = -pert(l) (g2) + -npcrt(l) ( g 2) 

An 



A) 



1 A 2 
ln(Q 2 /A 2 ) + A 2 -Q 



(3) 



^ote that V2C 2 = C 2 CSS [1]. 
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Figure 1. (a) Sudakov form factor versus transverse separation b for three Q 2 values: 
Qi = 2 GeV, Q 2 = 5 GeV, and Q 3 = 10 GeV, with all & = 1/2, and where we have set C\ = 2e~~< E , 
C 2 = e" 1 / 2 and A QCD = 0.242 GeV. The dotted curve shows the result obtained with a^ s , and 
Aqcd = 0.2 GeV for Q 2 = 5 GeV, using the same set of Cj. In that case, evolution is limited by the 
(renormalization) scale [iR = t = {mscx^/xy Q, C\/b}, as proposed in [10], albeit the enhancement 
at small 6-values due to the quark anomalous dimension is not neglected here, (b) Saturation 
behavior of pion's electromagnetic form factor, calculated in the AFS at NLO with commensurate 
scale setting (see text) and including a mass term (with m q = 0.33 GeV) in the BHL ansatz [11] 
for the soft pion wave function. The scheme parameters are defined in (17). Here b cr denotes the 
integration cutoff over transverse distances in (23). The momentum transfer values are as in part 
(a) 



where here and below A = Aqcd is the QCD scale parameter. 

To leading order in analytic perturbation theory (APT), one has 

T H (x, y, b, Q; /i R ) = 8C F a s an (/4)^o {^ybQ) , 

where C F = (N 2 - 1)/2N C = 4/3 for SU(3) C . The amplitude 



(4) 



V n (x, b, P ~ Q, Ci, C 2 , jj) = exp 



-s (x, b, Q, Ci, C 2 ) - s (x, b, Q, C u C 2 ) 



P w (x,6,C7i/6) 



(5) 



describes the distribution of longitudinal momentum fractions of the qq pair, taking into 
account the intrinsic transverse size of the pion state and comprising corrections due to soft 
real and virtual gluons, including also evolution from the initial amplitude V n (x, b, C\jb) 
at scale C\jb to the renormalization scale fi oc Q (more details and references are relegated 
to [9]). The main effect of the absence of a Landau pole in the running coupling a^ n is 
to make the functions s (x, b, Q, C±, C 2 ), s (x, b, Q, C±, C 2 ) well-defined (analytic) in the IR 
region and to slow down evolution by extending soft-gluon cancellation down to the scale 
Ci/b ~ Aqcd, where the full Sudakov form factor acquires a finite value, modulo its Q 2 
dependence (see LHS of Fig. 1). In addition, as we shall see below, the Sudakov exponent 
contains power-behaved corrections in (Ci/6A) 2p and (C 2 /t;QA) 2p , starting with p — 1. Such 
contributions are the footprints of soft gluon emission at the kinematic boundaries to the 
non-perturbative QCD regime, characterized by the transversal (or IR) and the longitudinal 
(or collinear) cutoffs. 

The pion distribution amplitude evaluated at the (low) factorization scale C\jb is ap- 
proximately given by 
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Table 1. Values of parameters entering the pion wave function [9]. The values in parentheses 
refer to the case m q = and the subscript "as" on /3 2 to the asymptotic distribution amplitude 



Input parameters Determined parameters 

m q = 0.33 GeV A = \- 10 - 01 (g ' 6 ) 

U = 0.1307 GeV 0^ = 0.871 GeV" 2 (0.743 GeV -2 ) 

(k 2 ) 1/2 = 0.352 GeV (0.367 GeV) 
Pqq = 0.306 (0.250) 



V n (x,b,C 1 /b,m q ) ~ 



2JW r 



0^ (x,Ci/b) E (x,b,m q ) . 



(6) 



To model the intrinsic transverse momenta of the pion bound state, we have to make 
an ansatz for their distribution. (For a recent derivation from an instanton-based model, 
see [12]). Here, I employ the Brodsky-Huang-Lepage (BHL) ansatz [11] and parameterize 
the distribution E(x, kj_,m q ) in the intrinsic transverse momentum k± (or equivalently the 
intrinsic inter-quark transverse distance b) in the form of a non- factorizing in the variables 
x and k± (or x and b) Gaussian function: 



V n (x, k_L,Ci/b, m q 



UV 2 

2N r 



$ w (x,Ci/6)E(x, k ± ,m q ), 



where 



$^ (x,d/b) = A$ as (x) = A6x(l-x) 



(7) 



(8) 



is the asymptotic distribution amplitude, with A being an appropriate normalization factor, 
and where 



with 



and 



E (x, k ± , m q ) = 167r /3 w g(x)T, (x, k ± ) E (x, m q ) 



E(x,k ± )=exp -^kig(x) 



E (x, m q ) = exp -(3 w m q g(x) 



(9) 



(10) 



(11) 



By inputting f n and the value of the quark mass m q and using g(x) = 1/ (xx), with x = 
(1 — x), we determine the parameters (we refer for more details to [9]) A, (3%, P qq , and 
(k^) 1/2 , tabulated in Table 1. 

We have now to calculate the Sudakov contribution within the AFS. Generically, the 
Sudakov form factor Fs (£, b, Q, C±, C2), i.e., the exponential factor in front of the wave 
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function, will be expressed as the expectation value of an open Wilson (world) line along a 
contour of finite extent, C, which follows the bent quark line in the hard-scattering process 
from the segment with direction (four-momentum) P to that with direction P' after being 
abruptly derailed by the hard interaction which creates a "cusp" in C. It is to be evaluated 
within the range of momenta from C\jb (IR cutoff) to C 2 ^Q (longitudinal cutoff) (where 
£ = x, x, y, y) and the region of hard interaction of the Wilson line with the off-shell photon 
is factorized out. Then the Sudakov functions, entering (5), can be expressed in terms of 
the momentum-dependent cusp anomalous dimension of the bent contour [1,13-15] to read 



1 r c ^Q da 

s (£, b, Q, &, c 2 ) = - / Jt r cusp ( 7 , <C(n)) (12) 

£ JCi lb u 



1 r c ^Q dfi 

'd/b \X 

with the anomalous dimension of the cusp given by 



r cusp ( 7 ,cCM) =2\n(^^ ] jA(aT^))+B(a: i \ f ,)) , 

ppert I pnpert /i q\ 

— 1 cusp ' 1 cusp ' V 1,3 / 

7 = ln(C2^Q/fj) being the cusp angle, i.e., the emission angle of a soft gluon and the bent 
eikonalized quark line after the external (large) momentum Q has been injected at the cusp 
point by the off-mass-shell photon, and where in the second line of (13) the superscripts 
relate to the origin of the corresponding terms in the running coupling. The functions A 
and B are known at two-loop order: 



7^(ar(/i))+/3(^)^/C(C 1 ,a s an (^)) 



and 



CF <^)) + W )c /<toM)V (14) 

7T 2 y 7T J 



B «») = -\ [K (d, < n (//)) + Q (e, C 2 , < n (/i))] 



The first term in (14) is universal, 2 while the second one as well as the contribution termed 
B are scheme dependent. The K-factor in the MS scheme to two- loop order is given by 
[1,13,16] 

K (Ci) = - C A - ^n f T F + A, In (16) 



2 In works quoted above, the cusp anomalous dimension is identified with the universal term, 
whereas the other (scheme and/or process dependent) terms are considered as additional anomalous 
dimensions. Here this distinction is irrelevant. 
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with Ca = Ac = 3, rif = 3, Tp = 1/2, and 7e being the Euler-Mascheroni constant. A set 
of constants Cj, (i = 1,2,3), which eliminate artifacts of dimensional regularization while 
practically preserving the matching between the re-summed and the fixed-order calculation, 
are [9] 

C 1 = 2exp(- 7E ), C 2 = exp(-l/2), C 3 = 2 exp (- 7e ) , C 4 = exp (-4/3) , 

AT = 4.565, k = 0. (17) 

The quantities /C, Q in (15) are calculable using the non-Abelian extension to QCD [1] 
of the Grammer-Yennie method for QED or employing the Wilson (world) lines ap- 
proach [13-15]. The soft (Sudakov-type) form factor depends only on the cusp angle which 
varies with the inter-quark transverse distance b ranging between C\/b and C 2 t;Q- The 
corresponding anomalous dimensions are inter-linked through the relation 2r cusp (a an (/i)) = 
7k « n (/i)) with r cusp (a s an (/i)) = C F a s an (/i 2 )/7r, which shows that \ lK = A(a s an (/i)). (Note 
that 7g = -7^ and 7q « n (At)) = -afV 2 )/ 71 "-) 

The leading contribution to the Sudakov functions s (£, b, Q, C±, C 2 ) (where £ = x,x, y, y) 
within our framework, is obtained by expanding the functions A and B in a power series 
in a an and collecting together all large logarithms ( ) In ( ^r^bQ) , which correspond to 



large logarithms In 
from the expression 



in transverse momentum space. The leading contribution results 



2 JCi/b pL I V /I 



a 



an(2) 



A w 



TV 



+ 1 



a 



an(l) 



7Y 



A® (d) 



(18) 



where the two-loop expression [8] for the strong coupling is to be used in front of A^\ 
whereas the other two terms are to be evaluated with the one- loop result. Let me remark at 
this point that in the following we ignore the difference between the analytic strong coupling 
squared and its "analytized" second power. These issues will be considered elsewhere. The 
specific values of the coefficients are 

A™ = C F , 
A^{C X ) = 1 -C F K(C 1 ) , 



2 / C 2 p 2 7e-i~ 



(19) 



in which the term proportional to A^ represents the universal part. The universal part of 
the Sudakov factor in LLA and including power corrections, reads 



jriuniv 



(a*f,Q) 




n in/i F /A 2 11 14 U \l4 / 



Li, ^ 



(20) 
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where Q represents the scale d(,Q and the IR matching (factorization) scale //p varies with 
the inverse transverse distance b, i.e., /ip = d/b- Note that the four last terms in this 
equation originate from the non-perturbative power correction (cf. (13)), and that Li 2 is the 
dilogarithm (Spence) function which comprises power-behaved corrections of the IR-cutoff 
(6A) and the longitudinal cutoff (Q/A). To complete the discussion about the Sudakov 
factor, I display the result obtained by neglecting power corrections: 



s (£,b,Q, d,d) 



Po 



(2A«g + J B«)ln|-2A« (Q-S) 



- -A™ 



x In 



Q Q-b s 



1 r 



+ 2 



In 2 (2Q) - In 2 (2b) 



|Aw{ln2-^[l + ln(2S) 



(21) 



where the convenient abbreviations [10] Q = In and b = In ^ have been used. Note 
that expressions given in the literature by other authors are erroneous. 

In the following, (18) is evaluated numerically to NLLA with appropriate kinematic 
bounds [9] to ensure proper factorization at the numerical level. The electromagnetic pion 
form factor in next-to-leading logarithmic order has the following form in LO of Tr: 



2 /"i fi r°° 

F n (Q 2 ) = -A 2 nC F f* dx dy bdba^ (/x R ) $ v {x)$M 
6 Jo Jo Jo 

b 2 (xx + yy) 



x exp 



cxp 



<rn\ - + - 
q V xx yy 



K (V^Qb) 



x exp [-S (x, y, b, Q, d, C 2 , d)} , 
whereas in NLO it reads 

f v (q 2 ) = i6A 2 ttC f (t^=) l ldx Jo dy l™ hdha T (/4) *M*M 

b 2 (xx + yy) 



(22) 



x exp 



cxp 



pW q i — + — 



1 



1 



n q 



K(^y-Qb) 



x exp (-5 (x, y, 6, Q, Ci, C 2 , C 4 )) 1 1 + [/ UV (x, y, Q 2 //4) 

+/ir (a;, V,Q 2 /l4) + fc (x, y )] | . 

In these equations the Sudakov form factor, including evolution, is given by 

S (x, y, b, Q, Ci, d, d) = s (x, b, Q, d, d) + s (x, b, Q, d, d) + (x ^ y) 

-8r(d/b,fi R ) 



(23) 



(24) 



with the "evolution time" [9] 
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Figure 2. Space-like pion form factor calculated within the AFS. Further details are provided 
in the text. LO calculation (dashed-dotted line); NLO calculation (dotted line). The dashed line 
gives the result for the soft, Feynman-type contribution, computed with m q = 0.33 GeV in the 
pion wave function, and the solid line represents the sum of the NLO hard contribution and the 
soft one [9]. The data are taken from [19,20]. 



r 



Ci 



rv 2 dk 2 < n(1) (fe 2 ) 
Jc'l/b 2 k 2 



An 



1 ^ lnQu 2 /A 2 ) 
A> n \n(C 2 /(bA) 2 



1 

Wo 



In 



(Ci/by 



-In 



I/. 2 - A 2 



£±. _ a 2 



(25) 



and the functions fa taken from [17]. I present predictions for F n in Fig. 2, adopting the 
BLM commensurate-scale method [18], and setting fj,p = C\/b and /xblm = AtRexp(— 5/6), 
where /x R = C 4 f(x, y)Q = C 4y /xyQ. 

As one sees, the hard contribution to F n (Q 2 ) within the AFS and with a BLM-optimized 
choice of scales provides a sizeable fraction of the magnitude of the form factor - especially 
at NLO. No artificial rising at low Q 2 of the hard form factor appears, as in conventional ap- 
proaches, so that this region is dominated by the Feynman-type contribution [21]. Moreover, 
the self-consistency of perturbation theory has been improved, as one infers from the satura- 
tion behavior of the scaled form factor, presented on the RHS of Fig. 1. Indeed, Q 2 F n (Q 2 ) 
accumulates the bulk of its magnitude below b CT AQcr>/Ci < 0.5, i.e., for short transverse 
distances, where the application of perturbative QCD is sound. Even better predictions can 
be obtained, using a more accurate pion distribution amplitude, recently derived in [22] with 
QCD sum rules and non-local condensates. 
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2.2 Power Corrections to Pion Form Factor 



The rationale of global analyticity, i.e., analyticity as a whole, implies 



Q 2 F n (Q- 



Cdx £dy $r (v, Q 2 ) Th(x, y, Q 2 , a s (Q 2 )W:(x, Q 2 ) 

A f dx f dyxyxy[T n (x,y,Q 2 ,a s (Q 2 ))] &n 
Jo Jo 



(26) 



wherre A is a normalization constant for the pion distribution amplitude, taken again to 
be the asymptotic one. Without the analytization requirement, the pion form factor is not 
Borel-summable (see, e.g., [23]), but only an asymptotic series in the coupling constant. 
Analytization entails 



< (Q 2 ) 



d£ 



it Jo £ + Q 2 — it 



(27) 



where the spectral density is the dispersive conjugate of all powers n of a s . For the 

leading-order expression under consideration the spectral density becomes [8,24,25] 



p (Q 2 ) = lma s (-Q 2 ) 



7T 



1 



Pi ln 2 (Q 2 /A 2 ) + 7T 2 



(28) 



Then (27) reduces to 



(V 



1 



dz 



2ixi Jc z — Q 2 + ie 



a s (z), 



(29) 



where C is a closed contour in the complex z-plane with a branch cut along the negative real 
axis, assuming exactly the form of (3), as proposed by Shirkov and Solovtsov [8]. Recasting 
the strong coupling in the form 



a. 



1 1 f°° / \ 



(30) 



with the plus sign corresponding to the case \z\/A 2 > 1 and the minus one to \z\/A 2 < 1, 
and inserting it into (26), we find after some standard manipulations the Borel transform of 
the scaled pion form factor at leading perturbative order [4] : 



Q 2 F n (Q 2 ) 



(i) 

an 



daexp (-<t/3i lnQVA 2 )^^) 



(i) 

an ' 



(31) 



The Borel image of the form factor reads 



7T JO Jo JO 



di 



C + xy 



x 



Q 2 



A 2 



IV 



-t 



(32) 
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This expression has no IR renormalons in contrast to approaches that use the conventional 
one-loop a s parameterization [23]. 

Hence, the integration over the Borel parameter a can be performed without any ambi- 
guity to arrive at the following result for the pion form factor 



Q 2 F„ (Q< 



(i) 

an 



16nC F A 



1 r 1 



Pi Jo 



dw 4>(w) 



+ 



i _ wQ 2 
1 A 2 



(33) 



The remaining integration can be carried out analytically to arrive at an expression derived 
in [4]. Here I only display the expression for the physically relevant case Q 2 3> A 2 : 



Q 2 F^ (Q 2 )]^ = IGttCWL \-a s (Q 2 ) + O (a 2 s ) 



1 A 2 



In 2 



A 2 



- 2 In 



(34) 



referring for further details to [4]. 



3. POWER CORRECTIONS TO DRELL-YAN PROCESS 

As a second example of the AFS, I discuss the derivation of power corrections to the 
inclusive Drell-Yan cross-section with the large scale Q 2 being here the invariant lepton pair 
mass. Citations to previous works and full details of the derivation are given in [4]. Consider 
the logarithmic derivative of the unrenormalized expression of the eikonalized Drell-Yan cross 
section, with the notations of [2] : 

rflnH/ DY ^ n(1) ^ 



rflnQ 2 



^V^^ 8 ^^'"" 1 )- (35) 



The following important remarks are now in place: (i) The argument of a s is taken to 
depend on k± to account for higher-order quantum corrections, originating from momentum 
scales larger than this [16]. (ii) The integral over the transverse momentum is not well- 
defined at very small mass scales owing to the Landau singularity of the QCD running 
coupling in that region, (iii) The evaluation at the edge of phase space is sensitive to the 
regularization applied to account for power corrections due to soft gluon emission transient 
to nonperturbative QCD. 

Imposing analytization as a whole and integrating over transverse momenta, we obtain 

[nW (Q 2 )] an = J™dae-^ l <^ 2 )flV(a) (36) 

with a Borel transform given by 
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n<2(<r)= 



4C F 



2z,2* 



7T 



sin (irafli 



,2 A 2 \ -2<r^i 



+ 



6 2 A 



'b 2 A 2 



6 2 A 2N 



where 



^(0 



d 2 " 2e g 1 1 



(e" 2 ** 6 - l) 



(37) 



(38) 



(2vr) 2 - 2e g 2 g 2 + £ 

Combining denominators in Eq. (38) and carrying out the integrations over £, we then find 



C, 



V /Jan 7T V 'JO 



-<t/3i ln(4/6 2 A 2 ) _ 



X 



X 



a Pi + e 



T(l-aPi-e) + J2 



T(l + aPi) 

;-i) n 



(b 2 A 2 \ n+1 ~ al3l ~ e 



n + 1 — <7/5i — e 



(39) 



with f{b 2 A 2 /4) being a complicated expression, provided in [4] and r(x, y), denoting the 
incomplete Gamma function. The first term in (39), viz., the integral over a, diverges 
for a Pi = 0, i.e., for small values of a s (k±) (or equivalently for large transverse momenta 
k±). This UV divergence is regulated dimensionally within the MS renormalization scheme 
adopted here. Were it not for the terms containing powers of bA, expression (39) and that 
found by Korchemsky and Sterman [2] (namely, their equation (18)) would be the same. 
In our case, however, the imposition of analytization cures all divergences related to IR 
renormalons that are generated by the T-functions whenever a Pi is an integer different from 
zero. 



Let us concentrate on the second term in 



IlW (Q 2 ) that gives rise to power corrections. 



Retaining only the leading contribution in b A , we find 

, 2a2 x 6 2 A 2 5 2 A 2 b 2 A 2 



O (6 4 A 4 ) 



(40) 



with the constant coefficients [4]: a = 0, ai ~ 3.18, a 2 ~ —2.51. Now one can expand 
the integral in the first term of IT^ (Q 2 ) in powers of 6 2 A 2 and regulate the UV pole at 

L J an 

a Pi = dimensionally. For o~Pi an integer, both terms inside the bracket have poles, but 
they mutually cancel so that their sum is singularity-free and the integral finite. Retaining 
terms of order 6 2 A 2 , the main contribution stems from the leading renormalon at a Pi = 1: 



[IK') = [IK') (<3 2 ) 

with the perturbative part being defined by 

[n« (q 



PT 



+ 



pow 



C F , ln (C/b 2 A 2 ^ 
In 



pt nPi ln (Q 2 /A 2 ) ' 



(41) 



(42) 
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a result coinciding with that obtained in [2]. Power corrections in the impact parameter b 
are encoded in the second contribution (b 2 A 2 Cl): 



n« (Q 2 )] pow = S + b 2 S 2 (b 2 A 2 ) + O (6 4 A 4 ) , 



(43) 



where 



and 



5 = — 5-ao = 



5 2 (6 2 A 2 ) = -%AVln^-a 2 

v ' Ati(3i 4 



(44) 



(45) 



The DY cross-section ITdy, comprising the leading logarithmic perturbative contribu- 
tion (Sudakov exponent SW) and the first power correction (in 6 2 A 2 ) reads (with the (in- 
dependence arising due to collinear interactions) 



where 



W BY (b, Q) = exp -S PT (b, Q) - b 2 S 2 (b, Q) + ... 



S 2 (b, Q) ~ S 2 (b 2 A 2 ) InQ + const. 



(46) 



(47) 



The Sudakov factor, representing the perturbative tail of the hadronic wave function, sup- 
presses constituent configurations which involve large impact space separations, while the 
exponentiated power corrections in b 2 , being of nonperturbative origin, provide enhancement 
of such configurations (since S 2 (b 2 A 2 ) is always negative). Hence, the net result is less sup- 
pression of the DY cross-section and also enhancement of the pion wave function in b space 
with the endpoint region 6A ~ 1 being less enhanced relative to small b transverse distances. 



4. CONCLUSIONS 

I have presented a theoretical framework, based on analytization that enables the calcu- 
lation of perturbative gluonic corrections (Sudakov form factor), as well as power-behaved 
ones that are linked to nonperturbative effects in QCD. Moreover, one can calculate the ab- 
solute normalization of the power corrections to hadronic observables systematically without 
any renormalon ambiguity from the outset. 
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